O 

(N 



(N 



X 



Direct-Photon Production in p+p Collisions at ^/s = 200 GeV at Midrapidity 

A. Adare," S. Afanasiev,^^ C. Aidala,^^ N.N. Ajitanand,^^ y. Akiba,47.48 h. Al-Bataineh ,^2 J. Alexander,^^ 
K. Aoki,30.47 L_ Aphecetche,55 J. Asai,'''^ E.T. Atomssa,^! R. Averbeck,^^ T.C. Awes,''^ ^ Aznioun,^ 
. Babintsev,2i M. Bai,^ G. Baksay/"^ L. Baksay/"^ A. Baldisseri," K.N. Barish/ P.D. Barnes,^^.^ B. Bassalleck,"*! 
A.T. Basye,^ S. Bathe/ S. Batsouli,''^ V. Baublis,''^ C. Baumann,^^ A. Bazilevsky,^ S. Belikov,6'EI 
R. Bennett,^'* A. Berdnikov,5° Y. Berdnikov,^'^ A.A. Bickley," J.G. Boissevain,33 H. Borel,^'^ K. Boyle.^"* 
M.L. Brooks,33 H. Buesching,^ V. Bumazhnov,^! G. Bunce,^'''^ S. Butsyk,^^ C.M. Camacho,^^ S. Cainpbell,^'* 
B.S. Chang,63 W.C. Chang,^ J.-L. Charvet," S. Chernichenko,^! C.Y. Chi/^ M. Chiu,22 I.J. Choi,63 
L.K. Choudhury,^ T. Chujo,^^ P. Chung,53 A. Churyn,^! V. Cianciolo,*^ Z. Citron;^" B.A. Cole/^ p. Constantin,33 
M. Csanad/*^ T. Cs6rg6,62 T. Dahms,'^'^ S. Dairaku,30>47 k. Das/** G. David/ A. Denisov/i D. d'Enterria/i 

A. Deshpande/s- E.J. Desmond/ O. Dietzsch/^ A. Dion/"^ M. Donadelli/i O. Drapier/^ A. Drees/'' 
K.A. Drees/ A.K. Dubey/^ A. Durum/i D. Dutta/ V. Dzhordzhadze/ Y.V. Efremenko/^ F. Ellinghaus/^ 

T. Engelniore/2 A. Enokizono/^ H. En'yo/^-''* S. Esumi/^ K.O. Eyser/ B. Fadem/^ D.E. Fields/^'*^ 
[. Finger/ M. Finger, Jr./ F. Fleuret/^ S.L. Fokin/^ Z. Fraenkel/^'R J.E. Frantz/'* A. Franz/ A.D. Frawley/^ 



T. Engelmore/^ A. Enokizono/^ H. En'yo/'-*" S. Esumi,='» K.O. Eyser/ B. Fadeni/» D.E. Fields/"' 
M. Finger/ M. Finger, Jr./ F. Fleuret,^^ S.L. Fokin,^^ Z. Fraenkel,^^'El J.E. Frantz,^'* A. Franz/ A.D. Frawley, 
K. Fujiwara,-*^ Y. Fukao/O'-*^ T. Fusayasu,"" I. Garishvili,^^ A. Glenn," H. Gong,^" M. Gonin/i J. Gosset,!" 
Y. Goto/^'-^s R. Granier de Cassagnac,^! N. Grau/^ S.V. Greene/^ M. Grosse Perdekanip,22.48 t. Gunji/o 
H.-A. Gustafsson/5'EI A. Hadj Henni/^ J.S. Haggerty,^ H. Hamagaki,i° R. Han,'*^ E.P. Hartouni,32 K. Haruna,20 
E. Haslum/5 R. Hayano/° X. He/^ M. Heffner/2 T.K. Hemmick,54 T. Hester,^ J.C. Hill,25 M. Hohlmann/^ 
W. Holzniann/3 K. Honima/° B. Hong,28 T. Horaguchi/O' ^T, 57 ^ Hornback,^^ g Huang/^ T. Ichihara,"*^' 
R. Ichimiya/'^ H. linuma/O' '''^ Y. Ikeda,^^ K. Imai,^^'^^ j. Imrek/^ M. Inaba/^ D. Isenhower/ M. Ishihara,"*^ 
T. Isobe/0'47 M. Issah/3 A. Isupov/^ D. Ivaniscliev,46 B.V. Jacak,54,[t| j_ jia/2 j^^i2 gj^^ Johnson,^ 
K.S. Joo,39 D. Jouan/'' F. Kajihara,!'' S. Kametani,''^ N. Kamihara,"* J. Kamin/'' J.H. Kang,^^ j_ Kapustinsky,33 
D. Kawall,36'48 a.V. Kazantsev/^ T. Kempel,25 A. Khanzadeev/^ K.M. Kijima/° J. Kikuchi,^^ B.I. Kim/^ 
D.H. Kim/9 j3 J Kim/3 g Kim/2 g jj Kim/3 ^ Kinney," K. Kiriluk/i A. Kiss/^ E. Kistenev/ J. Klay,32 
C. Klein-Boesing,37 L. Koclienda,''^ B. Komkov,''^ M. Konno,^* J. Koster,22 A. Kozlov,^! A. Kral,i3 a. Kravitz,i2 
G.J. Kunde,33 K. Kurita,47^49 Kurosawa,*^ M.J. Kweon,28 Y. Kwon.^^^ G.S. Kyle,42 j^^ Lacey,53 y.S. Lai,i2 
J.G. Lajoie,25 D. Layton,22 A. Lebedev,25 D.M. Lee,33 K.B. Lee,28 T. Lee,52 ^ j Leitcli,33 M.A.L. Leite,^! 
B. Lenzi,5i X. Li,^ P. Liebing,48 T. Liska,i3 A. Litvinenko/^ H. Liu,'*2 M.X. Liu,33 B. Love,^^ D. Lyncli,^ 

C. F. Maguire,^^ Y.I. Makdisi,^ A. Malakliov,26 M.D. Malik/i V.I. Manko,29 E. Mannel,i2 Y. Mao/^'^^ 
L. Masek,8.24 h. Masui,^^ F. Matathias,i2 M. McCumber,54 p.L. McGaugliey,33 N. Means,^'' B. Mereditli,22 

Y. Miake/* P. Mikes,24 K. Miki,^^ A. Milov,^ M. Mislira,3 J.T. Mitchell,^ A.K. Mohanty,-* Y. Morino,i° 
A. Morreale,"^ D.P. Morrison,*^ T.V. Moukhanova,^^ D. Mukliopadliyay,^^ J. Murata,'''^' ''^ S. Nagamiya,^'^ 
J.L. Nagle," M. Naglis,^! M.I. Nagy,!^ I. Nakagawa,''^- Y. Nakamiya,20 T. Nakamura,20 K. Nakano,47.57 
J. Newby,32 M. Nguyen,54 T. Niita,^^ R. Nouicer,^ A.S. Nyanin,29 E. O'Brien,^ S.X. Oda,!" C.A. Ogilvie,25 
M. Oka,5S K. Okada,^* Y. Onuki,47 A. Oskarsson,35 M. Oucliida,20 K. Ozawa/° R. Pak,^ A.P.T. Palounek,33 
V. Pantuev,23. 54 y. Papavassiliou,'*2 Park,52 ^ Park,28 S.F. Pate,'*2 H. Pei,25 J.-C. Peng,22 H. Pereira," 
V. Peresedov/*^ D.Yu. Peressounko,29 C. Pinkenburg,^ M.L. Purschke,^ A.K. Purwar,33 H. Qn,^^ J. Rak/^ 
A. Rakotozafindrabe,3i I. Ravinovich,^^! K.F. Read,'*3.56 g Rembeczki," K. Reygers,37 V. Riabov,'*^^ Y. Riabov,'*^ 

D. Roacli,59 G. Roclie,34 S.D. Rolnick,'^ M. Rosati,25 S.S.E. Rosendahl,35 P. Rosnet,34 P. Rukoyatkin,26 
P. Ruzicka,24 V.L. Rykov,^^ B. Salilmueller,37 N. Saito,30.47,48 gakaguclii,^ S. Sakai,^^ K. Sakashita,"*^' 

V. Samsonov,''^ T. Sato,^** S. Sawada,^^ K. Sedgwick,'^ J. Seele," R. Seidl,22 A.Yu. Semenov,25 V. Semenov,2i 
R. Scto,'^ D. Sharma/i I. Shein,2i T.-A. Shibata/^' K. Sliigaki,20 M. Shimomura,^^ K. Slioji,30.47 p_ ghukla,'' 
A. Sickles,'^ C.L. Silva,^i D. Silvermyr,43 C. Silvestre/"* K.S. Sim,28 B.K. Singli,3 CP. Singli,3 V. Singli,3 
M. Slunecka,^ A. Soldatov/i R.A. Soltz/2 W.E. Sondheim,33 S.P. Sorensen,^^ I.V. Sourikova,^ F. Staley," 
P.W. Stankus,43 ^_ Stenlund,35 M. Stepanov,42 A. Ster,62 g_p^ ^^.^^^6 Sugitate,20 C. Suire,44 A. Sukhanov/ 
J. Sziklai,'^2 E ]y[_ Takagui,^! A. Taketani,47. 48 Tanabe,^^ Y. Tanaka,''° K. Tanida/^- ^^2 ^.J. Tannenbaum,^ 
A. Taranenko/3 P. Tarjan,!^ H. Tliemann,54 T.L. Thomas/^ M. Togawa,30'47 a. Toia,54 L. Tomasek,24 
Y. Tomita,^^ H. Torii/O'47 R.g_ Towell,i V-N. Tram,3i I. Tserruya,^! Y. Tsucliimoto,20 C. Vale,25 H. Valle,^^ 
H.W. van Hecke,33 A. Veicht,22 J. Velkovska,^^ R. Vertesi,!^ A.A. Vinogradov/^ M. Virius,i3 V. Vrba/* 
E. Vznuzdaev/6 X.R. Wang,'*2 y. Watanabe,''^' "^^ F. Wei/^ J. Wessels,^^ S.N. Wliite,^ D. Winter,i2 C.L. Woody,^ 
M. Wysocki/i W. Xie,48 Y.L. Yamaguclii,^" K. Yamaura,2" R. Yang,22 A. Yanovicli,2i J. Ying,!^ S. Yokkaichi/^' ^8 
G.R. Young,43 I. Younus/i I.E. Yushmanov/^ W.A. Zajc,i2 Q. Zaudtke,37 C. Zliang,43 S. Zhou,^ and L. Zolin26 



X : T. Isobe/'0'47 

CD 
I 

Ph. 



> 

m 
cn 
in 
in 

vn 
o 

(N 



2 



(PHENIX Collaboration) 

'Abilene Christian University, Abilene, Texas 79699, USA 
^Institute of Physics, Academia Sinica, Taipei 11529, Taiwan 
^Department of Physics, Banaras Hindu University, Varanasi 221005, India 
^Bhabha Atomic Research Centre, Bombay 400 085, India 
^Collider-Accelerator Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
^Physics Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
''University of California - Riverside, Riverside, California 92521, USA 
^Charles University, Ovocny trh 5, Praha 1, 116 36, Prague, Czech Republic 
^Science and Technology on Nuclear Data Laboratory, China Institute of Atomic Energy, Beijing 102413, P. R. China 
Center for Nuclear Study, Graduate School of Science, University of Tokyo, 7-3-1 Hongo, Bunkyo, Tokyo 113-0033, Japan 

University of Colorado, Boulder, Colorado 80309, USA 
'^Columbia University, New York, New York 10027 and Nevis Laboratories, Irvington, New York 10533, USA 
Czech Technical University, Zikova 4, 166 36 Prague 6, Czech Republic 
'^Dapnia, CEA Saclay, F-91191, Gif-sur-Yvette, Prance 
'^Debrecen University, H-4010 Debrecen, Egyetem ter 1, Hungary 
'^ELTE, Eotvos Lordnd University, H - 1117 Budapest, Pdzmdny P. s. 1/A, Hungary 
'"'Florida Institute of Technology, Melbourne, Florida 32901, USA 
'^Florida State University, Tallahassee, Florida 32306, USA 
'^Georgia State University, Atlanta, Georgia 30303, USA 
^"Hiroshima University, Kagamiyama, Higashi-Hiroshima 739-8526, Japan 
^'IHEP Protvino, State Research Center of Russian Federation, Institute for High Energy Physics, Protvmo, 142281, Russia 
"'^^^ University of Illinois at Urbana- Champaign, Urbana, Illinois 61801, USA 
Institute for Nuclear Research of the Russian Academy of Sciences, prospekt 60-letiya Oktyabrya 7a, Moscow 117312, Russia 
^'^Institute of Physics, Academy of Sciences of the Czech Republic, Na Slovance 2, 182 21 Prague 8, Czech Republic 

^^lowa State University, Ames, Iowa 50011, USA 
Joint Institute for Nuclear Research, 141980 Dubna, Moscow Region, Russia 
^^KEK, High Energy Accelerator Research Organization, Tsukuba, Ibaraki 305-0801, Japan 
^^Korea University, Seoul, 136-701, Korea 
Russian Research Center "Kurchatov Institute", Moscow, 123098 Russia 
^"Kyoto University, Kyoto 606-8502, Japan 
^'Laboratoire Leprince-Ringuet, Ecole Poly technique, CNRS-IN2P3, Route de Saclay, F-91128, Palaiseau, France 
Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
^'Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 
^^LPC, Universite Blaise Pascal, CNRS-IN2P3, Clermont-Fd, 63177 Aubiere Cedex, France 

Department of Physics, Lund University, Box 118, SE-221 00 Lund, Sweden 
Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA 
^^Institut filr Kemphysik, University of Muenster, D-48149 Muenster, Germany 
^^Muhlenberg College, Allentown, Pennsylvania 18104-5586, USA 
Myongji University, Yongin, Kyonggido 449-728, Korea 
^''Nagasaki Institute of Applied Science, Nagasaki-shi, Nagasaki 851-0193, Japan 
University of New Mexico, Albuquerque, New Mexico 87131, USA 
^^New Mexico State University, Las Cruces, New Mexico 88003, USA 
'^^Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 
^^IPN-Orsay, Universite Paris Sud, CNRS-IN2P3, BPl, F-91406, Orsay, France 
Peking University, Beijing 100871, P. R. China 
^''PNPI, Petersburg Nuclear Physics Institute, Gatchina, Leningrad region, 188300, Russia 
^''RIKEN Nishina Center for Accelerator-Based Science, Wako, Saitama 351-0198, Japan 
^^RIKEN BNL Research Center, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
'^^ Physics Department, Rikkyo University, 3-34-1 Nishi-Ikebukuro, Toshima, Tokyo 171-8501, Japan 
Saint Petersburg State Polytechnic University, St. Petersburg, 195251 Russia 
^' Universidade de Sao Paulo, Instituto de Ftsica, Caixa Postal 66318, Sao Paulo CEP05315-970, Brazil 

^^Seoul National University, Seoul, Korea 
^^Chemistry Department, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 
^'^ Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, New York 11794-3400, USA 
^^SUBATECH (Ecole des Mines de Nantes, CNRS-IN2P3, Universite de Nantes) BP 20722 - 44307, Nantes, Prance 

'''"University of Tennessee, Knoxville, Tennessee 37996, USA 
^''Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152-8551, Japan 
^^Institute of Physics, University of Tsukuba, Tsukuba, Ibaraki 305, Japan 
^^Vanderbilt University, Nashville, Tennessee 37235, USA 
^"Waseda University, Advanced Research Institute for Science and 
Engineering, 17 Kikui-cho, Shinjuku-ku, Tokyo 162-0044, Japan 
"' Weizmann Institute, Rehovot 76100, Israel 



3 



Institute for Particle and Nuclear Physics, Wigner Research Centre for Physics, Hungarian 
Academy of Sciences (Wigner RCP, RMKI) H-1525 Budapest II4, POBox 49, Budapest, Hungary 
^■'Yonsei University, IPAP, Seoul 120-749, Korea 
(Dated: March 4, 2013) 

The differential cross section for the production of direct photons in p+p collisions at ^/s = 200 
GeV at midrapidity was measured in the PHENIX detector at the Relativistic Heavy Ion Collider. 
Inclusive-direct photons were measured in the transverse momentum range from 5.5-25 GeV/c, 
extending the range beyond previous measurements. Event structure was studied with an isolation 
criterion. Next-to-leading-order perturbative-quantum-chromodynamics calculations give a good 
description of the spectrum. When the cross section is expressed versus xt, the PHENIX data 
are seen to be in agreement with measurements from other experiments at different center-of-mass 
energies. 

PACS numbers: 25.75.Dw 



I. INTRODUCTION 

Direct photons are defined as photons that do not orig- 
inate from hadronic decays. In hadron-hadron collisions, 
direct photons at large transverse momentum {pt) are 
predominantly produced by the fundamental quantum- 
chromodynamics (QCD) 2-to-2 hard-scattering subpro- 
cesses, g + (7 — !■ 7 + g and q + q — > 7 + 5, where the former 
subproccss, which dominates in p+p and A-l-A collisions, 
is called "the inverse QCD Compton effect" [1]. This sub- 
process is one of the most important of the QCD 2-to-2 
subprocesses for three reasons: 

1. the photon is the only outgoing particle in funda- 
mental QCD 2-to-2 subprocesses that is a single 
particle, which can be measured to high precision; 

2. the scattered quark has equal and opposite trans- 
verse momentum to the direct-photon so that the 
transverse momentum of the jet from the frag- 
mented quark is also precisely known (modulo kx 
or multisoft gluon effects 0); and 

3. it is directly sensitive to the gluon distribution func- 
tion of the proton times the distribution function 
of quarks, which is precisely measured in deeply 
inelastic lepton-proton scattering. 

If both the direct photon with px and rapidity and 
the away side jet at yj are detected then, to the extent 
that the 9 + g — > 7 + g subprocess can be neglected in 
p+p collisions due to the predominance of gluons over 
anti-quarks, the jet opposite to the direct photon is a 
quark Q . Tagging jets with direct photons provides an 
excellent method of studying any medium effect on the 
energy or fragmentation of the outgoing quark. Further- 
more, the cross section for 5 + (7 7 -|- (7 in LO pQCD [[1] 
in scattering of hadron A from hadron B takes on the 



simple form for the reaction A + B^j + q + X: 

, 2 1 ^ J = xigA{xi) F2b{x2) 
dp^ dyj dyj 

naasiQ^) fl + cos9* 2 



(1) 



1 + cos ( 



3s2 V 2 

+X2gB{x2) F2a{xi) 
iraasiQ'^) ( I- cos 9* 2 



3s2 V 2 ' l-cosr 
where the parton kinematics are fully determined by 



Xl — XT ■ 



-y-i -i_ p~yj 



X2 = XT ■ 



, (2) 



and XT = 2ptI \fs. The parton-parton cm. energy \fh — 
^x\X2S, where yfs is the A + B cm. energy; the energy 
of the direct photon in the parton-parton cm. system is 



P; ^E* = Vs/2, where 



Pt ^Pt 



■sine* 



(3) 



and cos 9* = tanh(?/.^ — yj)/2 is the cm. angle of the out- 
going 7 with respect to hadron A. In EqlU ^^(xijQ^) 
and gB{x2,Q^) are the gluon structure functions of 
hadron A and hadron B. At leading order F2a{xi,Q^) 
and F2b{x2, Q^) are structure functions measured in deep 
inelastic scattering (DIS) of e-|-^, given by F2k{x, Q^) = 
^Yla fai^j Q^); whcrc fa{x, Q^) are the distributions 
in the number of quarks of type a, with electric charge 
Ca (in units of the proton charge) in hadron A^ . In hard 
g + q — >■ 7 + 9 scattering in p+p collisions, the struck 
quark is 8 times more likely to be a u quark relative to a 
d quark. For production in nuclei, the ratio is somewhat 
less according to the ratio of the atomic number to the 
atomic mass. 



^ The relation F2 = F2{DIS) is true only in leading order pQCD. 
'Deceased They can be different in higher order pQCD. But those difference 

t PHENIX Spokesperson: | j acak® skipper . physics . sunysb. edu| is accounted in the theory. 
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Beyond leading order, direct photons can be produced 
either by bremsstrahlung from any quark hne in a 2-to- 
2 subprocess, e.g. g + q^g + q + ^, or in a parton 
shower from fragmentation that forms a jet. In both 
these cases the photon is accompanied by jet fragments 
so that observing photons isolated from jets enhances the 
contribution from the fundamental 2-to-2 subprocesses. 
Naturally, all these effects must be taken into account 
in theoretical calculations of direct photon production in 
pQCD, and such calculations [1] are generally in excel- 
lent agreement with all previous measurements, includ- 
ing those from PHENIX However, decreasing the 
uncertainties of both measurement and theory and ex- 
tending the range to larger pr is desirable. Measure- 
ments with and without an isolation criterion allow more 
specific comparisons of theoretical models and a better 
understanding of photons coming from bremsstrahlung 
and parton fragmentation. 

Measurement of direct photons in p-\-p collisions pro- 
vides an important baseline for measurements in heavy- 
ion collisions. Once produced, a photon emerges from the 
reaction almost unaffected since it only interacts electro- 
magnetically. Initial state modifications of the distribu- 
tion functions in nuclei can be accessed by measurements 
in p+A or d+A collisions. Similarly, direct photons pro- 
vide a reference free from final state effects on the out- 
going quark, which at LO initially balances the pT of the 
direct photon. 

In this paper, we report a major update of the direct 
photon cross section measurement in p+p collisions at 
y/s = 200 GeV. The present data has more than an or- 
der of magnitude improved statistics than that reported 
in 0. It has extended the highest pT reach of the mea- 
surement from 15 GeV/c to 25 GeV/c. We compare the 
data to pQCD calculations and other direct photon data 
in hadronic collisions. 

The paper is organized as follows: Sec. HI] describes 
the experimental setup. Sec. IIIII describes the analysis 
method. Results are in Sec. IIVI followed by a discus- 
sion (Sec. |V| and summary (Sec. |VI|. The measured 
invariant cross sections are tabulated in the Appendix. 



II. EXPERIMENTAL SETUP 

Photons were detected in the PHENIX central arm 
detectors by two electromagnetic calorimeter (EMCal) 
arms (West, East) each covering tt/2 rad in azimuthal 
angle (0) and \r]\ < 0.35 in pseudorapidity. Each arm 
is divided into 4 sectors in azimuth. All 4 sectors in 
the West arm are lead scintillator sampling detectors 
(PbSc). In the East arm, 2 sectors are PbSc and 2 
sectors are lead glass Cerenkov detectors (PbGl). The 
sectors are composed of independent towers with granu- 
larities of At? X = 0.011 X 0.011 and 0.008 x 0.008 
for the PbSc and the PbGl, respectively. A tower con- 
tains ^ 80% of the photon energy hitting the center 
of the tower. During data taking, the relative gain of 



the detectors was monitored using a light pulser calibra- 
tion system. The absolute energy calibration was based 
on the known minimum-ionizing energy peak of charged 
tracks, energy-momentum matching of identified electron 
tracks and the measured value of the tt*^ — >■ 77 invari- 
ant mass. The linearity of the energy response was ob- 
tained from beam tests and the dependence of the mea- 
sured tt'^ mass on its momentum. The energy resolution 
was determined using the width of the ir'^ peak, and was 
as/E = 8.1%/v^(GeV) ® 5.7%. The systematic error 
on the absolute energy scale is less than 1.5%. The time 
of flight (ToF) as measured by the EMCal with a reso- 
lution of better than 1 ns was used to reduce the cosmic 
ray background. 

The dynamic range of the electronics is saturated for 
the highest energy clusters (~ 25 GeV) measured in this 
analysis. The size of this effect is estimated using a con- 
volution of the maximum energy limit of each tower (26 
GeV typical) and the fraction of energy deposited in the 
central tower of a electromagnetic shower cluster. This 
effect was found negligible compared to the statistical 
uncertainty for the very high energy photons that are 
affected. 

The drift chambers (DC) and the innermost layer of 
the pad chambers (PCI) provide charged track informa- 
tion and were used to veto charged hadron clusters in 
the EMCal. Hits in the beam-beam counters (BBCs) po- 
sitioned at pseudorapidities 3.1 < \r]\ < 3.9 were used to 
measure a collision vertex from the time difference be- 
tween hits in both BBCs, and to monitor the luminosity. 
The PHENIX detector is described in detail elsewhere 

a. 



III. ANALYSIS 
A. Event selection 

The results in this paper are based on the data sam- 
ple taken during the 2006 RHIC run. A high pT photon 
sample was collected with an EMCal trigger in which the 
analog sum of signals from a 4 x 4 adjacent set of EMCal 
towers was greater than a nominal energy threshold of 2 
GeV, which was in coincidence with the minimum-bias 
trigger, corresponding to ^8.0 pb^^ integrated luminos- 
ity. The integrated luminosity was determined from the 
rate of a minimum-bias trigger that required hits in the 
BBCs and a collision vertex within 30 cm of the nomi- 
nal center of the interation region. At -y/s = 200 GeV 
this trigger selects 23.0 nib of the inelastic p + p cross 
section. This was measured with 9.7% uncertainty us- 
ing a Vernier-scan technique 0, [1] ■ This corresponds to 
about 55% of the inelastic p+p cross section. For the 
BBC trigger rate of 250 kHz typical of the sample used 
in this analysis in the crossing rate of 9.4 MHz, the effect 
of multiple collisions per bunch crossing was < 2% and 
is neglected. 
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B. Photon selection 

Photon candidates were reconstructed from EMCal 
clusters within an EMCal fiducial volume defined to ex- 
clude the edge areas Ayy^O-lO and A0=O.1O wide, result- 
ing in a fiducial area of \rif\ < 0.25 in pseudorapidity and 
A(j)f = 7r/2 — 0.2 in azimuth for each of the two arms. 
Areas outside the fiducial volume were included when 
searching for a partner candidate for 7r° — >■ 77 decays to 
suppress n'^ background, and to measure activity around 
direct photon candidates (an isolation cut, section UlI Ep 
to suppress bremsstrahlung and fragmentation photons 
in the reconstructed photon sample. 

Photon candidates were required to have px > 5 
GeV/c and a ToF measured by the EMCal to be within 
±5 ns from the expected arrival time for the photons 
originating at the vertex. This requirement reduced the 
background from cosmic rays by an order of magnitude. 
The remaining contributions are estimated from the ToF 
distribution and corrected. The magnitude of this cor- 
rection was negligible in the region of px <~ 15 GeV/c 
and larger for higher px photons as the rate goes down 
(~ 8% contamination at px — 25 GeV/c). The back- 
ground from charged tracks was suppressed by requiring 
that the EMCal cluster shape be consistent with a sin- 
gle electromagnetic shower and that no charged track 
points to to the cluster. The shower shape cut efficiency 
for photons, evaluated using reconstructed tt*^ decay pho- 
tons, was 0.98 over the relevant px range. Most of the 
conversion e'^e~ pairs in the ~ 10% of radiation length 
of material between the DC and the EMCal are recon- 
structed as single photons because of the minimal mag- 
netic field in this region. Of course, no charged track 
would point to the cluster since the conversion would 
have happened after the DC. An additional 1% loss was 
attributed to these photon conversions from a GEANT |i9i| 
simulation with a reasonable input px distribution. The 
contribution of the other hadronic background (neutrons, 
K]^, albedo from magnet poles and other material) was 
studied with a detailed GEANT Monte Carlo simulation 
and found to be less than 1% of the photon sample. 

The fine granularity of the PHENIX EMCal resolves 
the two photons from tt'^ — >■ 77 decays up to tt^ px of 12 
GeV/c (17 GeV/c) in the PbSc (PbGl). A 50 % merg- 
ing probability corresponds to 7r° px of 17 GeV/c (25 
GeV/c). In the px range presented in this paper (up to 
25 GeV/c), merged photons can be separated from single 
photon showers in the EMCal and rejected using shower 
shape measurements with an efficiency > 90%. 



C. Direct photon signal extraction 



evaluated by a Tr'^-tagging method ^. In this approach 
the direct photon candidate was paired with each of the 
other photons in an event (a partner photon) to calculate 
the two-photon invariant mass M-^-y , which was required 
to be in the range from 105 to 165 MeV corresponding to 
±3(7 around the tt" peak (see Fig. [T]). Both photons were 
required to have a minimum energy -Emin = 0.5 GeV. 
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FIG. 1: Two-photon invariant mass distribution in the West 
arm where one of the photons has 5 < pt < 5.5 GeV/c. 

The combinatorial background under the tt'^ peak was 
evaluated and then subtracted by fitting the two-photon 
invariant mass distribution outside the peak region. The 
corrections for the underestimation due to photon con- 
versions and TT^ Dalitz decays were applied as a part of 
the partner photon efficiency. To avoid acceptance losses 
for tt" reconstruction, the edge areas of the EMCal out- 
side the fiducial region were included for partner photon 
selection as long as the primary photon was in the fidu- 
cial region. The minimum photon energy cut, the EM- 
Cal geometry, and inactive areas led to an underestimate 
of the 7r° decay photon yields. A correction for this was 
calculated using a single particle Monte Carlo simulation, 
which included the EMCal geometry, the configuration of 
dead areas, resolution, and the tt*^ spectrum shape from 
earlier measurements Figure [5] shows a px dependent 
multiplicative correction, denoted as 1 + R, to the tagged 
TT*^ photon sample. 1-l-R drops with increasing px when 
going from 5 GeV/c to 15 GeV/c due to the decreasing 
influence of the i?min cut. As the px increases further the 
correction stops decreasing due to an increasing merg- 
ing probability. The correction for asymmetric decays is 
more affected by the i?min cut, since symmetric decays 
start merging and are rejected from the photon sample 



In the obtained photon sample the majority of the 
background for direct photon measurements comes from 
decays of hadrons, primarily 7r° — 77 80%) and 
?7 — > 77 (~ 15%). The contribution from tt*' decays was 



^ In the previous measurement Q, we introduced a Tr'^-tagging 
method and a cocktail subtraction method. Both use statistical 
subtraction, but the Tr'^-tagging method uses the photon px on 
which it is easier to apply an isolation cut. 
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by the shower shape cut. Completely merged photons 
from high px tt'^ decays may resemble single photons in 
the EMCal; this residual contribution is corrected later. 
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FIG. 2: The correction factor to be applied to the number 
of tagged photons for the total contribution from vr" ((f+R) 
in Eq|4]) from a single particle Monte Carlo simulation. The 
error bar shows the systematic uncertainty. 



The contribution to the photon sample from hadronic 
decays other than from 7r°s was estimated relative to the 
7r° decay contribution based on the "q/n^ [l^, and uj/n^ 
and rj' /iT^ ratios from our measurements assuming 
rriT scaling (tested in [ll|). We denote the ratio of pho- 
tons from these decays to the photons from 7r° decays 
as A. At lower pt {pt ~ 5 GeV/c) this ratio, A, has 
a weak px dependence and approaches 0.235 as pr in- 
creases. Since the contribution of photons from tt*' to 
the background decreases for pT > 10 GeV/c due to 7r° 
photons merging, the value, A, starts to rise linearly at 
around px = 12 GeV/c, and at the highest px point at 
25 GeV/c it is 1.4 (0.94) for West (East) arm. 

The yield of direct photons iVdir was obtained from the 
inclusive photon yield iVjnci as follows: 
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FIG. 3: (color online) In the West arm, (a) number of photon 
clusters and (b) Background-photon fraction. From bottom 
to top, tagged tv° (N^o), total photons from unmerged tt" 
{{l + R)N^o), aU hadronic decay {{1 + A){1 + R)N^o), and the 
total background including an estimate of completely merged 
clusters. In the highest 3 bins where no Tr^-tagged photons 
were found, only the uncertainties of TV^q are shown. 
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Ndir = Nincl - (1 + A){1 + R)N^o, 



(4) 



FIG. 4: (color online) The same as Fig. [3] but for the East 



where N^o is the contribution from tt'^s evaluated with a 
tagging process. In this notation, {l + R)-N^o represents 
the total contribution from unmerged tt" — 77 decays, 
and A-{1+R)-Nt^(> is the contribution from other hadronic 
decays. 

Figures [3] and m show different contributions to the in- 
clusive photon spectrum separately for West and East 
spectrometers. In the highest px bins, where no tt'^- 
tagged photons were found (due to low statistics and high 
merging probability), N^o was set to 0]to to safely cover 
the other hadronic-decay channels. 



D. Direct Photon Cross Section 

Based on the extracted direct photon yields A^dir in 
each Apx wide bin in transverse momentum, the invari- 
ant cross section of direct photon production was calcu- 
lated as follows: 



(Pa _ 1 1 iVdir 1 1 

dp^ C 2Trpx Apx Ay e euc 
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where £ is the integrated luminosity for the analyzed 
data sample, Ay is the rapidity range, e corrects for the 
acceptance including the trigger live area, photon recon- 
struction efficiency, trigger efficiency and pt smearing 
due to EMCal energy resolution, ebias corrects for the fi- 
nite efficiency of BBCs to trigger on high pT events. The 
latter was measured from the ratio of t:^ yields obtained 
using the EMCal based high pT photon trigger with and 
without BBC trigger requirements, and was found to be 
0.78 ± 0.02, independent of pr- 

A single particle Monte Carlo simulation, which in- 
cluded the configuration of detector active areas and res- 
olutions, was used to evaluate the corrections for the ac- 
ceptance and Pt smearing. The small differences between 
the distribution of azimuthal position of reconstructed 
photons between data and simulation (Fig. [5] and [6]) 
served to estimate the systematic uncertainty. The ef- 
fects of Pt smearing were determined by varying the in- 
put Pt spectrum of the simulation. The same simulation 
framework was used to evaluate the propagation of the 
1.5% scale uncertainty in the EMCal energy measure- 
ments to the final direct photon spectrum. 

The terms from Eq. |3]used to calculate A^dir contribute 
to its uncertainty as follows (as tabulated in Table IVT)) : 



5N^ 



(6) 



where W is defined as A^inciZ-Adir- At low pT with high 
backgrounds (W » 1), the contributing uncertainties 
on A'dir are amplified by a factor W ot W ~1. In this pT 
range the dominant uncertainty comes from the correc- 
tion for untagged photons ((1 -I- R)) from 7r°-decay due 
to its sensitivity to the minimal energy cut, i^min- In the 
higher pT bins the biggest systematic uncertainty of the 
component is in the merging effect (which is included in 
the (1+R) term), since most of tt'^s are merged in the EM- 
Cal in this Pt region. However this effect on the direct 
photon is small, suppressed by a factor (Vl^ — 1) << 1, 
because of the small background fraction and approaches 
zero at highest pt- 



E. Study of the effect of isolation cut 

In this section we investigate the effect of an isola- 
tion cut on direct photons to determine the fraction 
coming from Conipton and annihilation processes, which 
are expected to be isolated from jet activity. The en- 
ergy around the photon candidate in a cone of radius 
r — y/ {S-qY + {5(j)Y = 0.5 was required to be less than 
10% of the photon energy, in order to pass the isolation 
cut. The cone size is determined by the hadron correla- 
tion in a jet measurement (e.g. Fig. 6 in 0). The total 
energy in the cone was constructed by summing the en- 
ergy of electromagnetic clusters in the EMCal and the 
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FIG. 5: (color online) Number of photons (5 <pt < 10 GeV/c) 
as a function of the azimuthal angle (West arm). The his- 
togram shows the MC result normalized by the total count. 
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FIG. 6: (color online) Number of photons (5 <pt < 10 GeV/c) 
as a function of the azimuthal angle (East arm). The lower 
(upper) half corresponds to PbGl (PbSc) sectors. The his- 
togram shows the MC result normalized by the total count 
(PbSc and PbGl sectors are done separately.) 



momentum of charged tracks in the tracking system. To 
be counted as part of the cone energy, the minimum EM- 
Cal cluster energy was set at 0.15 GeV and the minimum 
track momentum at 0.2 GeV/c, close to the lower limit 
for charged particle reconstruction in PHENIX. To avoid 
inclusion of misidentified tracks, which, due to decays or 
photon conversions may mimic high pT tracks, the maxi- 
mum momentum for the tracks to participate in the cone 
energy calculation was set to 15 GeV/c. 

Aside from direct photons, the isolated photon sample 
(^inci) includes background photons from 7r° and other 
hadron decays. For tt'^ decays we consider photons that 
have a partner photon reconstructed in the EMCal ac- 
ceptance (n^o), and those which satisfy the isolation cri- 
teria if the partner photon is masked out {N^°). More 7r° 
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photons pass the isolation cut without the partner pho- 
ton energy, so n^^S is a subgroup of N'^q. The latter was 
used to estimate the isolated photons from tt" with miss- 
ing partner due to the energy threshold ii^min or EMCal 
masked areas, by multiplying by the same missing photon 
fraction R introduced in Eq. |4l 

To estimate the contribution of other hadron decays, 
the isolated photon candidate from 7r° {N!^°) is scaled. 
With this procedure, the isolation cut efficiency from jet 
fragments is taken into account, however there is an ad- 
ditional rejection due to its own partner photon. A single 
particle MC for rjs was used to include this effect. In the 
case of Tjs, for the lowest pr sample, the partner photon 
can be out of the EMCal acceptance because of a large 
opening angle, thereby reducing the rejection power. As 
it goes to high pr, the rejection power due to the partner 
energy becomes constant. 

Similar to EqHl the isolated direct photon yield (N^^.) 
was calculated using; 

iV^° = Nl^°i - + Nl%°R) - A''°{1 + R)N'^^? (7) 

In addition to A in Eq|31 includes the isolation cut 
effect of the hadron's own partner photon as described 
above. Different contributors to the isolated photon sam- 
ple N-^°i are shown in Fig. [7] and [3 

Uncertainties are propagated according to: 
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{1 + A'''°){1 + R)Nl%° 

dir 
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and 



{{1 + A^'°){1 + R) - l)Nll° 



dir 



(8) 



Smoothed functions from fits to the data of the Wo, 
Wi, W2, and W3 parameters were used. The overall trend 
is the same as in the case of inclusive photon measure- 
ment. However at low pT the systematic uncertainties are 
smaller due to smaller contribution from hadronic decay 
photons in the isolated photon sample. Appendix Table 
I VIII summarizes the systematic uncertainties for isolated 
direct photon measurements. 
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FIG. 7: (color online) Components of isolated photon clusters 
(West arm). Solid line (with error): isolated direct photon 
signal (Nf^), Dashed line : photons from tt" with missing 
partner (N'^oR), Dotted line : Photons tagged as (n'^S)- 
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FIG. 8: (color online) Components of isolated photon clusters 
(East arm) 



To estimate the remaining systematic uncertainty, each 
component contributing to the yield measurements was 
moved up and down by ±1(T of its systematic uncertainty. 
Figure El shows the variation due to this change. The val- 
ues are tabulated in Appendix Table [V] The total sys- 
tematic uncertainty was taken as the quadratic sum of 
these variations. 



F. Isolation over inclusive-direct photon ratio 

By taking the ratio of isolated direct photons to the 
inclusive-direct photons, some uncertainties such as pho- 
ton efficiency and the luminosity measurement cancel. 



To evaluate the rejection power of the isolation cut 
on jet fragmentation, the same ratio was calculated for 
photons from tt" (= N^^° /N^o). Here the systematic un- 
certainty is from the correction of combinatorial back- 
ground. 
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FIG. 9: (color online) Effects in the ratio of isolated-direct 
photons for (black points) —la and (red points) +la of sys- 
tematic uncertainty change: (closed circles) vr" counts, (open 
squares) untagged vr" corrections, (closed squares) nonphoton 
contributions, (open circles) hadronic components other than 

^0 



IV. RESULTS 

Figure [TU] shows the inclusive- direct photon cross sec- 
tion where the bin width correction in the data are ap- 
plied to the vertical direction. The data are compared to 



using the CTEQ 6M 
and the BFGII parton 



a NLO pQCD calculation [1J-| 
parton distribution functions 
to photon fragmentation function (FF) [l9| for three dif- 
ferent renormalization and factorization scales, from bot- 
tom to top, /x=2pT, PT and pt/2. The deviation of the 
data from the calculation is shown in the bottom panel. 
The present data are consistent with that we previously 
published in The highest px reach of the data is ex- 
panded from 15 GeV/c to 25 GeV/c. For the pt range 
(8 < pt < 25 GeV/c) a power law fit to the data gives 
n = 7.08±0.09{stat)±0.1{syst) with x^/NDF = 8.3/10. 
In the fit, all systematic uncertainties are treated as cor- 
related. 

To demonstrate the purity of the signal as a function 
of photon Pt, ratios to tt" spectrum are taken. The mea- 
sured TT^ cross section [20| is parametrized by the form 
[pb] = 1.777 X 10^^ p-^-^^ as show in Fig. [Ill The 
systematic uncertainty shown with a band does not in- 
clude the overall normalization uncertainty. Figure [12] 
shows the ratio of both the direct photon signal and pho- 
tons from tt'^ to the fit. The dotted line in the figure is 
at 2/(8.22 — 1) = 0.277, which is the analytic expecta- 
tion for the ratio of tt*^ decay photons to tt'~' in case of a 
pure power law behavior of the tt*^. The fraction of the 
direct photon contribution gets higher as the transverse 
momentum increases. The contribution of photons from 
TT^s deviates from the analytic expectation at higher px 
because most of the merged clusters are rejected from the 
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FIG. 10: (a) Inclusive direct photon spectra compared with 
an NLO pQCD calculation for three different renormalization 
and factorization scales as mentioned in the text (The same 
calculation was used in [^). The error bars denote point- 
by-point uncertainties, the error brackets show pr correlated 
uncertainties, (b) Comparison of the data and the pQCD 
calculation. 



photon sample. The data are systematically lower than 
the analytical line even at the lowest px- But as shown 
in Figim the systematic uncertainty of the fit constant 
is on the order of 10%. If this uncertainty is taken into 
account, the agreement is reasonable. 

Figure [T51 shows the ratio of isolated-direct photons to 
inclusive-direct photons, and isolated over inclusive pho- 
tons from tt'^s. Since the background subtractions are 
done independently, this ratio can exceed unity due to 
the uncertainty. The data are compared to two theoreti- 
cal calculations, which include the same isolation criteria. 
The cone size of the isolation cut is larger than the area 
within the PHENIX central arm acceptance around the 
isolated photon candidate and leads to an underestimate 
of the energy in the cone. This effect was not corrected 
for; instead the same acceptance was included in the the- 
ory calculation. However the theory calculation assumes 
no dead areas in the acceptance. The theory calculation 
varies at most by 2% (90% 92%), when the effect of 
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FIG. 11: (color online) The ratio of tt*^ spectrum [S^l to the 
power law fit for the range 5 < pr < 20 GeV/c. The sys- 
tematic uncertainty shown by the box does not include the 
overall normalization uncertainty. 
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FIG. 12: (color online) The filled circles (7dir/7r°): the ra- 
tio of the direct photon cross section to the n'^ cross section 
fit. The crosses (7^0 /vr''): the ratio of the background photon 
from 7r° decay to the tt'^ cross section fit. The tt" cross section 
used for the denominator is the fit curve explained in the text 
and shown in Fig. [11] The dotted line shows the analytical 
expectation of 7.^0 /tt". See the text for details. The sys- 
tematic uncertainty on the tt" photon contribution includes 
acceptance, smearing, 7r° photon counting, and untagged tt" 
photon probability. (These uncertainties enter as A2, C, and 
D in Table EB. 



generally they agree within the systematic uncertainty. 
Explanations of possible discrepancy were thought to be 
underlying event activity as well as the contribution of 
photons from quark fragmentation, which are not con- 
sidered in the theory calculation. However a study with 
an event generator (pythia tune A 22]) did not show 
any drop in the low pT region for the direct photons, 
while the level of isolated photons from tt" decays was 
well reproduced 



V. DISCUSSION 

Figure [TU] shows good agreement between the data and 
NLO pQCD calculations. While the calculations at low 
Pt seem low, the correlated systematic uncertainties in 
the data are such that the difference is not significant. 

Figure [14] compiles this data and other measurements 
of direct photons in p+p or p + p collisions from both 
collider and fixed target experiments, over a broad range 
of collision energy. Note that most of the collider data 
except for that of PHENIX apply an isolation cut in their 
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FIG. 13: The effect of isolation cut on direct photons and 
decay photons. Solid circles: Ratio of isolated direct photons 
to inclusive-direct photons. The statistical uncertainties are 
shown as black error bars and the systematic uncertainties 
are plotted as shaded bars. The solid and dashed curves are 
NLO pQCD calculations with three theory scales for BFGII 
[19I ] and one scale for the GRV [2l[ parton to photon frag- 
mentation functions. Open circles: Ratio of isolated photons 
from Tv" decays to all photons from vr" decays. 



the dead area is included. 

The isolation cut causes a large suppression of photons 
from tt''. This is expected as the is accompanied with 
other fragmentation products. At high pr (jpr >^ 10 
GeV/c), the ratio of isolated direct photon to inclusive- 
direct photons is typically more than 90% and matches 
the expectation from NLO pQCD calculations. At low 
Pt, the data are below the theory calculation, although 



^ In the previous paper Q , it was claimed that the effect of under- 
lying events was large. At that time, the Monte-Carlo calculation 
was done only for direct photon process and the ratio was scaled 
up to the NLO calculation. For the present work, a full mixture 
of processes was generated, so the result can be directly com- 
pared with data, assuming that pythia reproduces the physics 
correctly. 
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FIG. 14: (color online) Various direct photon cross section measurements in p+p and p + p collisions scaled by ' vs 
XT = 2pt/\/s- The legend shows the experiment and the center of mass energy [GeV] in parenthesis. Table |l] shows the 
references. 



photon selection as listed in Table HI Also note that the 
PHENIX measurement [23] having the lowest pT points 
uses the virtual-photon method. The cross sections are 
shown as a function of of xt — '^Pr/y/s and scaled by 
the empirical value of -^/s""" with riee = 4.5 (Eq. [9]). 
The effective power, UeS, is primarily sensitive to the 



quantum exchange governing the reaction but also has 
sensitivity to scale-breaking. For measurements of single 
particle or sing le jet inclusive px distributions, this xt 
scaling [l^, |43| | provides a data driven test of whether 
pQCD or some other underlying subprocess is at work, as 
well as providing a compact quantitative way to describe 
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TABLE I: Experimental data for Fig. 1141 Points consistent with are excluded from the plot. Reference [2^ is a good review. 



System 


Experiment 




Et range 


ri or xf range 


iso cut 


Data 




[Ref.] 


[GeV] 


[GeV] 








points 


P+P 


CMS [24] 


7000 


22 - 210 


\v 


< 1.45 


yes 


11 


P+P 


ATLAS [25] 


7000 


50 - 300 


\v 


< 0.6 


yes 


8 


PP 


DO [26] 


1800 


10.5 - 108.4 


\v 


< 0.9 


yes 


23 


PP 


CDF [27, 28] 


1800 


12.3 - 114.7 


Iv 


< 0.9 


yes 


16 


P+P 


PHENIX [29] 


200 


1.7-3.3 


\v 


< 0.35 


no 


4 


P+P 


PHENIX (This report) 


200 


5.3-25 


\v 


< 0.25 


no 


18 


P+P 


PHENIX [5] 


200 


3.75 - 15 


\v 


< 0.25 


no 


16 


PP 


UAl [30] 


630 


17-90 


\v 


< 0.8 


yes 


16 


PP 


UAl [30] 


546 


17-46 


k/ 


< 0.8 


yes 


6 


PP 


UA2 [31J 


630 


15.9-82.3 


\v 


< 0.76 


yes 


13 


P+P 


RUO [32] 


63 


4.7-8.7 


k/ 


< 0.8 


yes 


7 


P+P 


R806 [33] 


63 


3.75 - 11.50 


k/ 


< 0.2 


yes 


14 


P+P 


R807 [34] 


63 


4.75 - 10.36 


\v 


< 0.7 


yes 


11 


P+P 


R108 [35] 


62.4 


5.37- 12.44 


\v 


< 0.45 


yes 


8 


P+P 


E706 [36] 


38.8 


3.8-11 




1. < 77 < 0.5 


no 


9 


P+P 




ox.u 


^ S — Q 

O.O iy 


-0.75 <ri< 0.75 




e 



P+P 


E704 [37] 


19.4 


2.6-3.6 


|a::F| < 0.15 


yes 


5 


P+P 


NA24 [38] 


23.8 


3.3-6 


-0.65 < 77 < 0.52 


no 


5 


P+P 


WA70 [39] 


23.0 


4.1-5.7 


\xf\ < 0.05 


no 


5 


P+P 


UA6 [40] 


24.3 


4.2-6.3 


-0.1 < 77 < 0.9 


no 


9 


PP 


UA6 [41] 


24.3 


4.2-5.7 


-0.1 < r; < 0.9 


no 


6 



the data using the effeetive index, Ucs^xt, \/s) 



E- 



dp^ PTdpTdyd(j) 



-, i=:G{xt) 



(9) 



where Ed'^a /dp'^ is the invariant cross section for inclu- 
sive particle production with transverse momentum pT 
at cm. energy ^/s and xt = '^PTl\fs- It is important to 
emphasize that the effective power, ncs{xT,\/s), is dif- 
ferent from the power n of the invariant cross section at 
any given value of -y/s. 

For pure vector gluon exchange, or for QCD without 
evolution of as and the structure and fragmentation func- 
tions, neff = 4 as in Rutherford scattering. However, due 
to the scale breaking in QCD, the measured value of n^s 
depends on the Xt value and the range of y/s used in the 
computation [4^. 

For inclusive-direct photon production in p+p colli- 
sions at midrapidity, if we assume xi = X2 = xt and 
{cos9*) = 0, then from Eq[Tl the xt scaled inclusive 
cross section in pQCD is approximately: 



v^"="- E 



Hnott-4) XrgpjxT, Q"^) F2p{xT, Q^)as(Q^) 

xrgpixr) F2p{xt) 



(10) 



where xxgp is the gluon momentum distribution function 
in the proton and i^2p is the proton structure function 
measured in DIS and we assume that the empirical value 
rioff — 4 = 0.5 takes account of the scale breaking effects. 
The XT scaling of all the available data, with some excep- 
tion at low y^, is impressive. As one goes to higher xt, 
the power of the invariant cross section becomes softer. 
Figure [TJ] gives the same information as the agreement 
with the pQCD calculations Q but in addition shows the 
validity of pQCD directly from the data by a simple but 
powerful scaling rule. 



VI. SUMMARY AND CONCLUSIONS 

The invariant differential cross section for the produc- 
tion of direct photons in p+p collisions at y/s = 200 GeV 
at midrapidity was measured. It extends the pt reach 
up to 25 GeV/c. An NLO pQCD calculation agrees well 
with the measurement, supporting the validity of such 
calculations. 

The effect of an isolation cut on the direct-photon 
cross section was measured to be negligible (<10%) in 



agreement with NLO theoretical calculations. The iso- 
lation cut enhances the 5 + 9 ~> 7 + 9 contribution and 
suppresses a possible background of single photons from 
bremsstrahlung or jet fragmentation. The main utility 
of the isolation cut is that it reduces the background of 
photons from hadronic decays by a significant factor of 
~ 60%. Furthermore, the data are an important refer- 
ence for interpreting direct-photon spectra in heavy-ion 
collisions. 
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Tables of the measured invariant differential cross sec- 
tion, the ratio of isolated to inclusive-direct photon, and 
systematic uncertainties. 
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TABLE II: Cross section of midrapidity inclusive-direct pho- 
ton production in p+p collisions at ^/s = 200 GeV as a func- 
tion of transverse momentum (pr)- Asymmetric statistical 
uncertainties occur in pr bins with no tagged vr" counts. 



Pt 


E(fo/dp'^ 


Stat- 


Stat-f 


Syst 


[GeV/c] 


[ph ■ GeV-'^ ■ c^] 








5.25 


1.14e+03 


3.04e-|-01 


3.04e-h01 


4.78e-h02 


5.75 


6.13e-|-02 


1.92e-H01 


1.92e+01 


2.21e+02 


6.25 


3.48e-h02 


1.27e-h01 


1.27e-h01 


1.01e-H02 


6.75 


2.31e-F02 


8.50e-h00 


8.50e-h00 


6.24e-h01 


7.25 


1.36e-F02 


6.12e-h00 


6.12e-h00 


3.13e-H01 


7.75 


9.29e-h01 


4.41e-h00 


4.41e-f00 


1.95e-f01 


8.25 


6.70e-h01 


3.22e-|-00 


3.22e+00 


1.34e+01 


8.75 


4.83e-h01 


2.45e-H00 


2.45e-F00 


9.18e+00 


9.25 


3.21e-h01 


1.89e-|-00 


1.89e+00 


6.10e+00 


9.75 


2.04e-h01 


1.46e-h00 


1.46e-h00 


3.68e-h00 


11.00 


9.81e-h00 


4.23e-01 


4.23e-01 


1.67e-H00 


13.00 


2.97e-h00 


1.89e-01 


1.89e-01 


4.75e-01 


15.00 


1.06e-h00 


9.85e-02 


9.85e-02 


1.69e-01 


17.00 


3.38e-01 


5.51e-02 


5.51e-02 


5.42e-02 


19.00 


1.73e-01 


3.37e-02 


3.37e-02 


2.77e-02 


21.00 


8.82e-02 


2.06e-02 


2.01e-02 


1.50e-02 


23.00 


4.22e-02 


1.52e-02 


1.41e-02 


7.18e-03 


25.00 


2.87e-02 


1.41e-02 


9.07e-03 


4.30e-03 


TABLE III: Ratio of isolated/inclusive-direct photon (Fig. 
I13|l. Upper(-I-) and lower bounds(-) on systematics can be 
different, and are listed separately. 


< Pt > 


Ratio 


Stat 


Syst+ 


Syst- 


[GeV/c] 











0.151 
0.145 
0.130 
0.102 
0.103 
0.089 
0.067 
0.052 
0.056 
0.064 
0.052 
0.057 
0.046 
0.097 
0.071 
0.039 
0.100 
0.000 



TABLE IV: Ratio of isolated/inclusive photon from tt'' (Fig. 


< Pt > 


ratio 


Stat 


Syst+- 


[GeV/c] 








5.22 


0.450 


0.003 


0.018 


5.72 


0.443 


0.003 


0.018 


6.23 


0.433 


0.004 


0.017 


6.73 


0.419 


0.006 


0.017 


7.23 


0.418 


0.007 


0.017 


7.73 


0.409 


0.009 


0.016 


8.23 


0.422 


0.012 


0.017 


8.73 


0.403 


0.014 


0.016 


9.23 


0.392 


0.017 


0.016 


9. /3 


U.3(0 


0.U2U 


U.OiO 


10.79 


0.373 


0.015 


0.015 


12.82 


0.372 


0.029 


0.015 


14.84 


0.449 


0.054 


0.018 


16.86 


0.416 


0.083 


0.017 


18.88 


0.619 


0.171 


0.025 


20.89 


0.500 


0.354 


0.020 


TABLE V: Components of systematic uncertainties on 
isolated /inclusive-direct photon ratio (expressed as a percent- 
age of the center value). Note that the upper and lower vari- 
ations are treated separately. 


(pt) piOJagging 


piO_miss_ratio 


non_vertex 


other _hadron 


5.2 25.1 / -16.6 


20.0 / -14.3 


4.0 / -3.5 


7.7 / -5.9 


5.7 21.6 / -15.0 


16.4 / -12.4 


3.8 / -3.3 


7.1 / -5.7 


6.2 18.7 / -13.5 


10.8 / -8.8 


3.6 / -3.2 


6.9 / -5.7 


6.7 13.8 / -10.8 


7.6 / -6.6 


2.8 / -2.5 


5.2 / -4.5 


7.2 13.1 / -10.3 


5.7 / -5.1 


2.8 / -2.5 


5.3 / -4.5 


7.7 10.7 / -8.8 


4.5 / -4.1 


2.4 / -2.2 


4.5 / -4.0 


8.2 8.4 / -7.2 


3.4 / -3.2 


1.8 / -1.7 


3.4 / -3.1 


8.7 6.8 / -6.0 


2.7 / -2.6 


1.5 / -1.4 


2.9 / -2.6 


9.2 6.6 / -5.8 


2.6 / -2.5 


1.5 / -1.4 


3.0 / -2.7 


9.7 7.0 / -6.1 


2.6 / -2.5 


1.7 / -1.6 


3.4 / -3.1 


10.8 5.3 / -4.8 


2.1 / -2.0 


1.3 / -1.2 


4.0 / -3.6 


12.9 3.8 / -3.6 


1.9 / -1.8 


1.0 / -0.9 


4.6 / -4.1 


14.9 2.9 / -2.7 


1.8 / -1.7 


0.7 / -0.7 


4.0 / -3.6 


16.9 3.5 / -3.3 


3.0 / -2.9 


0.9 / -0.9 


10.4 / -8.2 


18.9 2.0 / -1.9 


2.6 / -2.5 


0.4 / -0.4 


7.3 / -5.9 


20.9 0.9 / -0.8 


2.8 / -2.7 


0.2 / -0.2 


3.1 / -2.8 


22.9 1.0 / -1.0 


7.2 / -6.3 


0.3 / -0.3 


6.5 / -5.7 


24.9 0.0 / 0.0 


0.0 / 0.0 


-0.0 / -0.0 


0.0 / 0.0 



5.23 


0.658 


0.014 


0.217 


5.73 


0.690 


0.017 


0.193 


6.23 


0.764 


0.022 


0.176 


6.73 


0.730 


0.021 


0.124 


7.23 


0.793 


0.027 


0.127 


7.73 


0.813 


0.029 


0.106 


8.23 


0.750 


0.028 


0.074 


8.74 


0.742 


0.029 


0.059 


9.24 


0.806 


0.035 


0.064 


9.74 


0.911 


0.047 


0.073 


10.83 


0.863 


0.027 


0.060 


12.85 


0.945 


0.040 


0.057 


14.87 


0.916 


0.055 


0.046 


16.89 


1.082 


0.118 


0.119 


18.90 


1.017 


0.119 


0.081 


20.91 


0.975 


0.096 


0.039 


22.92 


1.114 


0.193 


0.100 


24.92 


0.909 


0.091 


0.000 



15 



TABLE VI: Systematic uncertainties of inclusive-direct photon measurement. The percentage uncertainties of each component 
is shown: (A) the global factor (quadrature sum of components Al, A2, and A3), (B) inclusive photon counts, (C) vr" tagging, 
(D) the factor for the total tt" counts, and (E) the factor for the hadron contribution other than tt") and the contribution to 
the direct photon signal (A*l, B*W, C*(W-1), D*(W-1), E*(W-1)). W is the ratio of iVind/A^'dir- Individual components are 
(Al) Energy scale error transformed to the cross section, (A2) Acceptance and smearing, (A3) ctbbc and BBC trigger bias, 
(Bl) nonvertex, neutral hadron subtraction, (CI) combinatorial background subtraction, (C2) loss for conversion and Dalitz 
decay, (Dl) -Bmin calibration, (D2) input tt" spectrum in the MC, (D3) tt" merge model (correction for the complete merging), 
(D4) tt" merge model (cluster shape parametrization in the MC), (D5) Geometry and trigger mask, (El) Ratio of all hadronic 
decay to -k" contribution, (E2) Isolation with own decay partner, and (E3) tt" merge correction for other hadron contribution. 



Pt 


1/W 


Al 


A2 


A3 


A*l 


Bl 


B*W 


CI 


C2 


C*(W-1) 


Dl 


D2 


D3 


D4 


D5 


D*(W-1) 


El 


E2 


E*(W-1) 


total % 


5.25 


0.13 


10 


3 


10 


14.46 


1 


7.49 


3 


1 


20.53 


4 


2 








1 


29.74 


2 





12.98 


41.71 


5.75 


0.15 


10 


3 


10 


14.46 


1 


6.47 


3 


1 


17.31 


4 


2 








1 


25.08 


2 





10.94 


36.04 


6.25 


0.18 


10 


3 


10 


14.46 


1 


5.7 


3 


1 


14.86 


3 


2 








1 


17.58 


2 





9.4 


29.32 


6.75 


0.2 


10 


3 


10 


14.46 


1 


5.09 


3 


1 


12.93 


3 


2 








1 


15.3 


2 





8.18 


26.52 


7.25 


0.22 


10 


3 


10 


14.46 


1 


4.6 


3 


1 


11.38 


2 


2 








1 


10.79 


2 





7.2 


22.97 


7.75 


0.24 


10 


3 


10 


14.46 


1 


4.19 


3 


1 


10.1 


2 


2 








1 


9.58 


2 





6.39 


21.47 


8.25 


0.26 


10 


3 


10 


14.46 


1 


3.85 


3 


1 


9.02 


2 


2 








1 


8.56 


2 





5.71 


20.28 


8.75 


0.28 


10 


3 


10 


14.46 


1 


3.57 


3 


1 


8.11 


2 


2 








1 


7.7 


2 





5.13 


19.31 


9.25 


0.3 


10 


3 


10 


14.46 


1 


3.32 


3 


1 


7.33 


2 


2 








1 


6.95 


2 





4.63 


18.53 


9.75 


0.32 


10 


3 


10 


14.46 


1 


3.1 


3 


1 


6.64 


2 


2 








1 


6.3 


2 


1 


4.7 


18.01 


11 


0.38 


10 


3 


10 


14.46 


1 


2.67 


3 


1 


5.27 


2 


2 





1 


1 


5.27 


2 


2 


4.71 


17.14 


13 


0.46 


10 


3 


10 


14.46 


1 


2.18 


3 


1 


3.73 


2 


2 





3 


1 


5 


2 


3 


4.25 


16.45 


15 


0.54 


10 


3 


10 


14.46 


1 


1.84 


3 


1 


2.66 


2 


2 


1.7 


4 


1 


4.44 


2 


4 


3.76 


15.92 


17 


0.63 


10 


3 


10 


14.46 


1 


1.59 


3 


1 


1.88 


2 


2 


3.1 


6 


1 


4.4 


2 


6 


3.76 


15.77 


19 


0.71 


10 


3 


10 


14.46 


1 


1.41 


3 


1 


1.29 


2 


2 


6.7 


8 


1 


4.41 


2 


8 


3.35 


15.6 


21 


0.8 


10 


3 


10 


14.46 


1 


1.26 


3 


1 


0.82 


2 


2 


27.4 


10 


1 


7.56 


2 


10 


2.63 


16.59 


23 


0.88 


10 


3 


10 


14.46 


1 


1.14 


3 


1 


0.44 


2 


2 


56.3 


12 


1 


7.93 


2 


12 


1.67 


16.62 


25 


0.96 


10 


3 


10 


14.46 


1 


1.04 


3 


1 


0.12 


2 


2 


102.8 


14 


1 


3.99 


2 


13 


0.51 


15.04 



TABLE VII: Systematic uncertainties of isolated-direct photon measurement. The percentage uncertainties of each component 
are as given in the Table IVTl caption, except that the contributions to the direct photon signal are A*l, B*WO, C*W3, C*W2, 
D*W1, and E*W1. WO, Wl, W2, and W3 are defined in Eq|8] Individual components Al through E3 are also as given in the 
Table rvTl caption. 



Pt 


WO 


Wl 


W2 


W3 


Al 


A2 


A3 


A*l 


Bl 


B*WO 


CI 


C2 


C*W3 


C*W2 


Dl 


D2 


D3 


D4 


D5 


D*W1 


El 


E2 


E3 


E*W1 


Total % 


5.25 


3.69 


4.16 


0.1 


2.62 


10 


3 


10 


14.46 


1 


3.69 


3 


1 


8.3 


0.3 


4 


2 








1 


19.07 


2 


0. 


1 





8.33 


26.92 


5.75 


3.09 


3.27 


0.09 


1.98 


10 


3 


10 


14.46 


1 


3.09 


3 


1 


6.27 


0.3 


4 


2 








1 


14.99 


2 


0. 


1 





6.55 


22.92 


6.25 


2.65 


2.62 


0.09 


1.53 


10 


3 


10 


14.46 


1 


2.65 


3 


1 


4.83 


0.29 


3 


2 








1 


9.79 


2 


0. 


,3 





5.29 


19.06 


6.75 


2.32 


2.13 


0.09 


1.2 


10 


3 


10 


14.46 


1 


2.32 


3 


1 


3.8 


0.29 


3 


2 








1 


7.96 


2 


0. 


,3 





4.3 


17.63 


7.25 


2.08 


1.76 


0.09 


0.96 


10 


3 


10 


14.46 


1 


2.08 


3 


1 


3.04 


0.28 


2 


2 








1 


5.27 


2 


0, 


,4 





3.58 


16.23 


7.75 


1.89 


1.47 


0.09 


0.78 


10 


3 


10 


14.46 


1 


1.89 


3 


1 


2.48 


0.28 


2 


2 








1 


4.42 


2 


0. 


,5 





3.04 


15.73 


8.25 


1.75 


1.25 


0.09 


0.65 


10 


3 


10 


14.46 


1 


1.75 


3 


1 


2.06 


0.27 


2 


2 








1 


3.76 


2 


0. 


,6 





2.61 


15.41 


8.75 


1.64 


1.08 


0.08 


0.55 


10 


3 


10 


14.46 


1 


1.64 


3 


1 


1.73 


0.27 


2 


2 








1 


3.24 


2 


0. 


7 





2.29 


15.18 


9.25 


1.56 


0.94 


0.08 


0.47 


10 


3 


10 


14.46 


1 


1.56 


3 


1 


1.48 


0.26 


2 


2 








1 


2.82 


2 


0. 


7 





1.99 


15.02 


9.75 


1.49 


0.83 


0.08 


0.41 


10 


3 


10 


14.46 


1 


1.49 


3 


1 


1.28 


0.26 


2 


2 








1 


2.49 


2 


0. 


,9 





1.82 


14.91 


11 


1.37 


0.63 


0.08 


0.3 


10 


3 


10 


14.46 


1 


1.37 


3 


1 


0.95 


0.24 


2 


2 





1 


1 


2 


2 




1 


1 


1.55 


14.77 


13 


1.27 


0.46 


0.07 


0.21 


10 


3 


10 


14.46 


1 


1.27 


3 


1 


0.66 


0.23 


2 


2 





3 


1 


1.94 


2 


1. 


,5 


2 


1.46 


14.73 


15 


1.22 


0.36 


0.07 


0.16 


10 


3 


10 


14.46 


1 


1.22 


3 


1 


0.51 


0.21 


2 


2 


1.7 


4 


1 


1.89 


2 


2, 


,3 


3 


1.53 


14.72 


17 


1.2 


0.29 


0.06 


0.13 


10 


3 


10 


14.46 


1 


1.2 


3 


1 


0.41 


0.19 


2 


2 


3.1 


6 


1 


2.12 


2 


3, 


,2 


4 


1.58 


14.75 


19 


1.18 


0.22 


0.06 


0.1 


10 


3 


10 


14.46 


1 


1.18 


3 


1 


0.32 


0.18 


2 


2 


6.7 


8 


1 


2.43 


2 


4, 


,3 


6 


1.71 


14.81 


21 


1.15 


0.16 


0.05 


0.07 


10 


3 


10 


14.46 


1 


1.15 


3 


1 


0.23 


0.17 


2 


2 


27.4 


10 


1 


4.73 


2 


5. 


,4 


8 


1.59 


15.34 


23 


1.12 


0.1 


0.05 


0.04 


10 


3 


10 


14.46 


1 


1.12 


3 


1 


0.14 


0.15 


2 


2 


56.3 


12 


1 


5.85 


2 


6. 


,5 


9 


1.14 


15.68 


25 


1.09 


0.05 


0.05 


0.02 


10 


3 


10 


14.46 


1 


1.09 


3 


1 


0.07 


0.14 


2 


2 


102.8 


14 


1 


5.48 


2 


7. 


,5 


10 


0.67 


15.51 



16 



H. Fritzsch and P. Minkowski, Phys. Lett. B 69, 316 
(1977). [24 
A. Adare et al. (PHENIX Collaboration), Phys. Rev. D 

82, 072001 (2010). ' [25 
A. Adare et al. (PHENIX Collaboration), Phys. Rev. C 
80, 024908 (2009). [26 
P. Aurenche, , J. P. Guillet, E. Pilon, M. Werlen, and 

M. Fontannaz, Phys. Rev. D 73, 094007 (2006). [27 
S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 
Lett. 98, 012002 (2007). [28 
K. Adcox et al. (PHENIX Collaboration), Nucl. Instrum. 
Meth. A 499, 469 (2003). [29 
S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 
Lett. 91, 241803 (2003). [30 
A. Adare et al. (PHENIX Collaboration), Phys. Rev. D 
79, 012003 (2009). [31 
GEANT 3.2.1 Manual (1994), CERN W5013, URL 
http : //wwwasdoc .web . cern. ch/wwwasdoc/pdf dir/ge ant ■pd|c i2 
A. Adare et al. (PHENIX Collaboration), Phys. Rev.^D 

83, 032001 (2011). [33 
A. Adare et al. (PHENIX Collaboration), Phys. Rev. D [34 
83, 052004 (2011). 

L. E. Gordon and W. Vogelsang, Phys. Rev. D 48, 3136 [35 

(1993) . 

L. E. Gordon and W. Vogelsang, Phys. Rev. D 50, 1901 [36 

(1994) . 

P. Aurenche, A. Douiri, R. Baier, M. Fontannaz, and [37 
D. Schiff, Phys. Lett. B 140, 87 (1984). 

P. Aurenche, R. Baier, M. Fontannaz, and D. Schiff, [38 
Nucl. Phys. B 297, 661 (1988). 

H. Baer, J. Ohnemus, and J. F. Owens, Phys. Rev. D 42, [39 
61 (1990). 

H. Baer, J. Ohnemus, and J. F. Owens, Phys. Lett. B [40 
234, 127 (1990). 

J. Pumplin et al., JHEP 0207, 012 (2002). [41 
L. Bourhis et al., Eur. Phys. J. C2, 529 (1998). 
A. Adare et al. (PHENIX Collaboration), Phys. Rev. D [42 
76, 051106 (2007). 

M. Gluck, E. Reya, and A. Vogt, Phys. Rev. D 48, 116 [43 
(1993). 

T. Sjostrand (1994), cERN-TH-7112. 

W. Vogelsang and M. R. Whalley, J. Phys. G G23, Al 



(1997). 

V. Khachatryan et al., Phys. Rev. Lett. 106, 082001 
(2011). 

G. Aad et al. (ATLAS Collaboration), Phys. Lett. B 706, 
150 (2011). 

S. Abachi et al. (DO Collaboration), Phys. Rev. Lett. 77, 
5011 (1996). 

F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 73, 
2662 (1994). 

F. Abe et al. (CDF Collaboration), Phys. Rev. Lett. 74, 
1891 (E) (1995). 

A. Adare et al. (PHENIX Collaboration), Phys. Rev. C 
81, 034911 (2010). 

C. Albajar et al. (UAl Collaboration), Phys. Lett. B 209, 

385 (1988). 

J. Alitti et al. (UA2 Collaboration), Phys. Lett. B 288, 

386 (1992). 

A. L. S. Angelis et al. (CMOR Collaboration), Nucl. 
Phys. B 327, 541 (1989). 

E. Anassontzis et al., Z. Phys. C 13, 277 (1982). 
T. Akesson et al. (Axial Field Spectrometer Collabora- 
tion), Sov. J. Nucl. Phys. 51, 836 (1990). 
A. L. S. Angelis et al. (CCOR Collaboration), Phys. Lett. 
B 94, 106 (1980). 

L. Apanasevich et al. (FNAL E706 Collaboration), Phys. 
Rev. D 70, 092009 (2004). 

D. L. Adams et al. (E704 Collaboration), Phys. Lett. B 
345, 569 (1995). 

C. De Marzo et al. (NA24 Collaboration), Phys. Rev. D 
36, 8 (1987). 

M. Bonesini et al. (WA70 Collaboration), Z. Phys. C38, 
371 (1988). 

G. Ballocchi et al. (UA6 Collaboration), Phys. Lett. 
B436, 222 (1998). 

G. Sozzi et al. (UA6 Collaboration), Phys. Lett. B 317, 
243 (1993). 

R. Blankenbecler, S. J. Brodsky, and J. F. Gunion, Phys. 
Lett. B42, 461 (1972). 

R. F. Cahalan, K. A. Geer, J. B. Kogut, and L. Susskind, 
Phys. Rev. D 11, 1199 (1975). 



